Magnetic excitations in a bond-centered stripe phase: 
Spin waves far from the semi-classical limit 
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Using a spin-only model, we compute spin excitation spectra in a bond-centered stripe state with 
long-range magnetic order. We employ a bond operator formalism, which naturally captures both 
dimerization and broken spin symmetry in a unified framework. At low energies, the spin excitations 
resemble spin waves, but at higher energies they are very similar to spin-1 excitations of isolated 
spin ladders. Our theory does well describe neutron scattering data [J. M. Tranquada et at, Nature 
429, 534 (2004)] on La2-a;Baa;Cu04, pointing towards bond order in this material. 



For a number of cuprate superconductors it has been 
established, most notably by neutron scattering tech- 
niques [1-3], that states with incommensurate spin and 
charge correlations, commonly referred to as stripes, ap- 
pear over a significant range of the phase diagram. While 
in some materials these correlations remain dynamic 
[4, 5], in others they become static and apparently coexist 
with superconductivity at lowest temperatures. The role 
of these stripes for superconductivity has been discussed 
extensively [6-10], but is at present not fully understood. 

Recent experiments [11] have mapped out the spin ex- 
citations in an ordered stripe phase of La2-xBa2:Cu04 at 
a hole doping of x = ^. In this paper we will present a 
consistent theoretical description of these neutron scat- 
tering data for all energies. 

For cuprates with dopings near x — ^ the spatial pe- 
riod of the charge order is found to be four lattice spac- 
ings. Furthermore, the period of the spin modulation 
equals twice the charge modulation period, i.e., the or- 
dering wavevectors obey Kc = 2Ks. The microscopic 
structure of stripes has not yet been fully elucidated. 
Regarding the symmetry of the charge modulation w.r.t. 
reflection on a Cu-Cu bond axis one distinguishes bond- 
centered and site-centered stripes [10]. Theoretical pro- 
posals based on doping a paramagnetic Mott insulator 
predict the coexistence of superconductivity with bond 
order [4, 12], compatible with bond-centered stripes. On 
the other hand, theories starting from the ordered anti- 
ferromagnet and favoring site-centered stripes have also 
been put forward. For states with spin and charge stripe 
order it has been suggested that the magnetism in the 
hole-poor regions resembles that of the undoped antifer- 
romagnet, with local (tt, tt) ordering; then the stripes act 
as antiphase domain walls in the magnetic background. 

The neutron scattering spectra on La2_a;Baa;Cu04 [11] 
indicate a linearly dispersing mode near the ordering 
wavevector at low energies, not inconsistent with a spin- 
wave picture. However, at higher energies conventional 
spin waves (as have been employed in Refs. 13, 14 to 
stripe-ordered phases) do not appear to describe the 
data, which instead are closer to what one expects for 
gapped two-leg ladders. 
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FIG. 1: Sketch of the bond-centered stripe state considered 
here. We use a two-dimensional (2d) spin-only model, where 
the effect of the hole-rich stripes is to mediate a weak cou- 
pling, Js < 0, between hole-poor regions, accross the hole-rich 
ones. The hole-poor regions are spin ladders, characterized by 
couplings J, J'. The unit cell consists of 8 Cu sites, i.e., 4 spins 
in the effective Heisenberg model. 



The purpose of this paper is to point out that there 
exists a simple unified description of the spin dynamics 
at all energies. The crucial point is the coexistence of 
dimerization and weak magnetic order in bond-centered 
stripes. The spin ordering in the hole-poor region is far 
from the semi-classical limit (this is the limit described 
by conventional spin waves), but instead is better under- 
stood as being close to a quantum phase transition where 
magnetic order disappears. We capture this physics us- 
ing a spin-only model for the hole-poor regions, assum- 
ing that the hole-rich stripes provide a weak coupling 
between them. For simplicity, we focus on stripe states 
with charge-order period of four sites, where both the 
hole-rich and hole-poor regions have width two. Thus we 
consider a Heisenberg model with a geometry shown in 
Fig. 1, with couplings J, J', Jg for rungs, legs, and across 
the stripes, respectively. J and J' are antiferromagnetic, 
whereas J a is much weaker and ferromagnetic, mimicking 
the antiphase property of the stripe. This model has a 
quantum paramagnetic singlet phase and a magnetically 
ordered phase with ordering wavevector Kg — (37r/4, tt) 
in units of the Cu square lattice. The transition will oc- 
cur at a critical Jf which depends on J, J'. For | J^j of 
order J, J' one expects a well-ordered state with negligi- 
ble dimerization where excitations can be described by 
conventional spin waves, however, here we are interested 
in Ja closer to the critical value, where quantum effects 
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are strong and spin-wave theory needs to be modified. 

Bond operator formalism. A natural language for spin 
dimerization is the bond operator representation for spins 
[15]. For two spins 1/2 spins Si^ S2 one introduces 
bosonic operators for creation of a singlet |0) = -^(1 f J, 

) — I il)) and three triplet states tl^,a = x,y,z above 
a fictitious vacuum |0): 4|0) = --^{\ TT) - I ID), 
4l0) = 75(1 tt) + I ID), m = ^(1 Ti) + I it)). Then 
the following representation is exact: 5*^2 = f ± 
tl^s — icaff'ytjjt^) ■ The physical Hilbert space is defined 
by the constraint s^s + t]j,a = 1- For our system, Fig. 1, 
we will use bond operators for each rung of the ladders. 
The labels "1" ("2") will be assigned to the left (right) 
spin of each rung. 

The simplest singlet ground state corresponds to a 
product state of the form \4>q) = s- 10), where i labels 
rungs (dimers). Magnetic ordering on top of the dimer- 
izcd background corresponds to condensation of one of 
the triplets. A suitable product state is given by [16, 17] 



1^0) ~ n(4+Ae^''''''i)|0>. 



(1) 



where Q' = (tt, tt) is the ordering wavevector on the 
dimer lattice. A non-zero value of A corresponds to a 
state with magnetic order: it has a condensation ampli- 
tude A of tz triplets, and interpolates between the sin- 
glet product state (A = 0) and the classical Neel state 
(A = 1). The ordered state has 8 Cu sites per unit cell, 
which translates into two dimers per unit cell of the spin- 
only model, i.e., the dimers form a bipartite lattice with 
"checkerboard" order. The quantity A can be understood 
as a variational parameter: minimizing the energy of |0o) 
with the full Hamiltonian gives 



A^ = 



2J' 



2 J' -Js + J 



(2) 



for 2 J' — Jj, > J, and A = otherwise. This defines the 
location of the transition point at the mean-field level, 
we comment on corrections beyond mean-field below. 

To describe fiuctuations around the product state, 
1^0), we perform a basis transformation [17]: 
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(4 + Ae^Q'^'4)/yrTA^, 
(_Ae'Q'^'4+4)/v/rTA^, 



'■ix ' 



ly 



t. 



(3) 



The product state (1) takes the form |(/)o) = Hi ^1 

The Hamiltonian can now be re-written in this new 
basis. The lengthy expression, containing bilinear, cubic, 
and quartic boson terms, can be found in Refs. 17, 19. 
[Linear boson terms, which can appear in principle, ex- 
actly vanish for the correct choice of A (2)]. The sim- 
plest non-trivial approach is a harmonic approximation 



[17-19], similar to the Holstein-Primakoff approach to 
spin waves. It is obtained by resolving the constraint as 
s = = {1 — iUa)^^'^, expanding the root, and keep- 
ing only bilinear terms in the i bosons. The resulting 
Hamiltonian is found as: 



n 



EBqa 
9 



{ilJ-^a+^-^-)-i^) 



Here, tqa are the modified basis operators which have 
been Fourier transformed with respect to the dimer mo- 
mentum q. The functions ^q, Bq contain all information 
about the geometry of the system. For our model we have 



1 + (l-hA2)2 

+J Y^jT^ COS Qy- Y 5^ 
Js 1-A2 

Bqx = J COS qy - —Y^-j^ COS Q:,, 



A 



qz 



l-A^ 2A2 

Jz 7^ +K- 



1 + A2 ' (1 + A2) 



B, 



qz. 



Js 



with K 



Bqz = {J'cOSQy ^ COS 

4J' - 2J, and Aa 



.(1-A2)2 
^^^1 + A2)2 



A R 

iqy — -^qx, -^qy — -^qx* 

A Bogoliubov transformation of (4) leads to H = 

Zqa ^ciaTl^Tqa whcrc uj"^^ = A"^^ - B^^. In the param- 
agnetic state, A = 0, the excitations are threefold degen- 
erate and correspond to gapped triplet modes [15-20]. 
As has been shown in Rcf. 17, for A = 1 one recovers 
standard linear spin- wave theory, with tx,ty being the 
transverse spin- wave modes, and representing a lon- 
gitudinal mode (which is dispersionlcss and carries zero 
spectral weight at A = 1). Note that - although we have 
two dimers per magnetic unit cell - three magnetic modes 
arc obtained, which is due to the sublattice symmetry 
of the dimer lattice. The formalism provides a consis- 
tent description of both zero-temperature phases includ- 
ing the transition (albeit with mean-field critical expo- 
nents), and is able to describe Goldstone modes which 
arise on the ordered side close to the magnetic quantum 
phase transition. It can also be extended for external 
magnetic fields [17, 19], and has recently been employed 
to describe the quantum phase transitions of the coupled 
dimer system TlCuCls (Ref. 19). 

So far, we have restricted ourselves to a harmonic ap- 
proximation of independent bosons. Clearly, this approx- 
imation overestimates the tendency to magnetic order, 
and underestimates the energy of the excitations (as the 
hard-core repulsion is relaxed). For the paramagnetic 
phase, it is possbile to incorporate the repulsion, e.g., 
with the Brueckner method [21]; however, we do not 
know of a consistent way to extend such an approach to 
the ordered phase with a smooth connection between the 
two. Thus, we continue to work with the independent 
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FIG. 2: Dispersions ti;(k) of the magnetic modes, calculated 
for vertical bond-centered stripes with charge period four, see 
Fig. 1, using the harmonic approximation (see text). Solid: 
Doubly degenerate transverse mode. Dashed: Longitudinal 
mode. Parameter values are J — J' = 100 meV, Js = — 0.06J, 
and the J renormalization (see text) has been used [20, 22]. 
For these parameters the "condensation" amplitude is A = 

0. 27 and the staggered magnetization A/a — 0.30; where A — 

1, Ms — 0.61 would correspond to the spin-wave result in the 
semi-classical (Neel) limit of a square- lattice antiferromagnet. 



boson picture, but take into account interaction effects 
by using renormalized parameters. This idea has been 
studied for the single spin ladder in Ref. 20, where it was 
shown that a quantitatively accurate description of the 
single-particle dynamics is possible. The most important 
renormalization is found in the q = contribution to the 
excitation energy (the momentum-independent term in 
Aq), which is Jron = 1.77 J instead of J (for A = 0) - 
this will be used below. (Different renormalizations for 
the Aq and Bq cannot be easily taken in account, as wc 
have to require Goldstone's theorem to be fulfilled.) 

Results. We focus on parameters J — J' and small 
Js, and replace J by Jicn in Eq. (5) [22, 23]. Sample 
dispersion relations are plotted in Fig. 2; note that the 
momenta k in the Cu lattice and q in the dimer lattice are 
related by = qx/4^, ky = Qy. The doubly degenerate 
transverse Goldstone mode and the longitudinal mode 
are clearly distinguished; at higher energies, both modes 
resemble the dispersion of triplets in a two-leg ladder. 

Next we calculate the T — susceptibility as mea- 
sured in inelastic neutron scattering. A physical 
spin operator can create both single-particle and two- 
particle excitations out of the ground state, according to 



^a(k) 



[4(q) + ia(-q)](i 



)/2-ieai3jJ2, 



qi ^/3^ 



qi)t^(qi)(l + e''^^)/2 (where we have set s = 1). Thus, 
the magnetic response will consist of sharp peaks from 
elementary excitations and a two-particle continuum. In 
the following, we restrict our attention to the one-particle 
contributions, as the continuum will be hard to detect ex- 
perimentally. Averaging over the neutron spin polariza- 
tions, we obtain x"()i, iu) as sum of S peaks with weights 
determined by various matrix element terms. Impor- 
tantly, there appear factors of 1 1 — e''^^ |^ in the matrix el- 
ements for the one-particle processes, which enhance the 
response near k = (tt, tt). Also, the weight of the longitu- 
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FIG. 3: (color) Neutron scattering intensity, x"(k,ti'): The 
panels show fixed energy cuts as function of momentum in the 
magnetic Brillouin zone. Parameter values are as in Fig. 2. 
The S peaks have been replaced by Lorentzians with width 
r = 0.2J, and the responses of horizontal and vertical stripes 
have been added. The background arising from the two- 
particle continuum is not shown. The figure can be directly 
compared to Fig. 2 of Ref. 11. 

dinal mode is suppressed by a factor of [(1— A^)/(l-|-A^)]^. 

In order to compare with the experiment of Ref. 11 
we add the contributions from horizontal and vertical 
stripes, and plot the result as function of the external 
momentum at fixed energy. A sample result is shown in 
Fig. 3 where we have broadened the 5 peaks to account 
for the experimental resolution. Interesting features are 
observed: For small energies (e.g. 6 meV) the response 
consists of four peaks representing the four cones of spin- 
wave modes. With increasing energy the cones widen, 
however, the outer part becomes suppressed in inten- 
sity due to matrix element effects. Thus, the spectrum 
at 36 meV is dominated by the inner part of the spin- 
wave cones [close to (7r,7r)]. For even higher energies 
(where the inter-ladder coupling becomes unimportant), 
the modes gradually change their character - towards 
triplet modes of a dimerized ladder - and the scattering 
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intensity forms a diamond which moves outward with in- 
creasing energy. (The neglect of boson interactions in our 
approximation somewhat overestimates the mode disper- 
sion near the zone boimdary this docs not quahtativcly 
influence the result for x". It also has to be kept in mind 
that the precise form of the dispersion at high energies 
depends on microscopic details like the presence of ring- 
exchange etc.) Overall, one observes a good agreement 
with the experimental data of Tranquada et al. [11]. 

Conclusions. We have discussed the spin excitations 
in a bond-centered stripe phase, where strong dimeriza- 
tion coexists with weak coUinear magnetic order. Using 
a bond operator approach we have been able to consis- 
tently describe the modes for all energies: The transverse 
modes resemble semi-classical spin waves at low energies, 
but cross over into triplet modes of dimerized ladders at 
higher energies. The crossover energy, associated with a 
saddle point in the mode dispersion, is only a fraction of 
J for systems with weak antiferromagnetism being close 
to a magnetic quantum phase transition. 

Two ingredients are crucial for the "dual" character of 
the modes: (i) strong quantum fluctuations (dimeriza- 
tion), and (ii) the presence of two energy scales for the 
dispersion, namely a large J' and a small Jg. Interest- 
ingly, a plaqucttc-ordcred state with 2d charge modula- 
tion obeys only (i), hence there is no crossover in the 
character of the mode [18]. On the other hand, spin 
waves for stripes [13, 14] only capture (ii), resulting in a 
high-energy spectrum different from ours. 

The agreement of our description with the experiments 
[11, 24] suggests that the spin ordering in La2-a;Baj;Cu04 
is far from the semi-classical limit, and a spin- wave de- 
scription may only be used at lowest energies. Instead, 
strong quantum effects related to dimerization dominate 
the spectrum at higher energies. These results strongly 
point to the existence of bond order. A likely scenario is 
that cuprate stripes at .x = ^ arc bond-ccntorcd (Fig. 1) 
instead of site-centered, as the latter situation would im- 
ply hole-poor regions with a width of three sites, which 
then favor spin order closer to the semi-classical limit 
[13, 14]. [Of course, more complicated ordered states, 
e.g., with strong dimerization realized in a site-centered 
stripe structure, cannot be excluded.] 

At doping levels different from number of 

modifications can be considered: Even in the absence of 
a modulation in the site-charge density a bond-ordered 
(i.e. dimerized) state is still possible [4, 12], with a 
magnetic excitation spectrum very much similar to the 
higher-energy part of our results. For a; < | , the stripe 
distance will become larger upon decreasing doping [3], 
which on the one hand moves the ordering wavevector 
closer to {n,7r) - this presumably lowers the energy of 
the strong peak seen at (7r,7r) [13, 25]. On the other 
hand the hole-poor regions become wider, which sta- 
bilizes magnetic order, thus reducing dimerization and 
quantum effects. Therefore, for small doping one can ex- 



pect the semi-classical spin-wave approach [13, 14] to be 
more appropriate for the whole energy range. Finally, 
the present picture can in principle also apply to incom- 
mensurate modulation periods, as seen in recent STM 
experiments [26], e.g., if a superstructure like an array of 
period-4 and period-5 stripes is realized. 
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